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We derive an exact solution for the total kinetic energy of noninteracting spinless electrons at 
half-filling in two-dimensional bipartite lattices. We employ a conceptually novel approach that maps 
this problem exactly into a Feynman-Vdovichenko lattice walker. The problem is then reduced to 
the analytic study of the sum of magnetic phase factors on closed paths. We compare our results 
with the ones obtained through numerical calculations. 
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Introduction. — Non- interacting tight-binding electron 
models at half-filling in two-dimensional (2D) bipartite 
(e.g., square and hexagonal) lattices have recently re- 
ceived jienewed attention due Xo their role in condensed 
mattera and particle jphysicsa. For instance, several 
quantum field theoriescl arise in a natural way from 2D 
tight-binding lattice fermion problems at half filling Jar 
hexagonalEm (e.g., 2D graphited) and square latticescm. 
These quantum field theoriesa are important to the prob- 
lem of dynamical symmetry breaking, which plays a cen- 
tral role in many current areas of research; for instance, 
they provide a possible mechanism for generating t 
fermion mass spectrum in elementary particle physic 
Furthermore, noninteracting 2D tight-binding electrons 
in a perpendicular magnetic field have been the subject of 
intense study in areas of current interest like mesoscopic 
structures and the quantum Hall effect. More recently, 
the behavior of the kinetic energy of a 2D non-interacting 
electron gas under the influence of both a periodic poten- 
tial and a magnetic fieldJias been analyzed by Hasegawa 
et al.u and many othersLJtl. These results have beerure- 
lated to mean-field approaches to the t — J modclEl of 
high-temperature superconductors. It is of value to ob- 
tain analytical results for this important problem that 
has recently motivated many perturbative and numerical 
studies. 

The goal of this paper is to present exact results that 
relate the kinetic energy of the half-filled Fermi sea of 
tight-binding electrons with sum-over-paths on the lat- 
tice. The problem is then reduced to the study of phase 
factors on closed paths. From the computational point 
of view, we would like to present an alternative to the 
standard approaches. From a physical point of view, and 
following Feynman's program, we would like to reformu- 
late this quantum problem as an "average over histories" . 

Recently, the lattice path approach has also been ap- 



ims; the computation of: 
I, and (2) the supercon- 



plicd to two very diverse probl 
(1) equilibrium crystal shapes 
ducting transition temperature, in wire micronetworks 
and Josephson-junction arraysta 

Fermi-sea energy. — The kinetic energy of spinless elec- 
trons on a 2D lattice in a uniform magnetic field is de- 
scribed by the Hamiltonian H = — X)<«j> c t c i ex P[*<^ij]i 
where < ij > refers to nearest-neighbor sites and 0y = 

2tt J j A-dl in units of the flux quantum. Also, \j >= 
Cj\0 > defines a localized state centered at site j, and we 
will work on the {\j >} basis. The above Hamiltonian has 
a fractal quantum energy spectrum as a function of an 
applied magnetic field and also plays a central role in the 
physics of particles with fractional statistics (anyons)O. 

At half-filling, the Fermi-sea ground-state energy is 
the sum of the lowest N/2 eigenvalues, where N is the 
total number of sites. Since the energy spectra of bi- 
partite (e.g., square and hexagonal) lattices are sym- 
metric under {E} —* {— E}, we can write the Fermi- 
sea ground-state energy per site as Et = jj J2e<o & = 
— (z/2N) Tr \Hq/z\. Here, Hq is the corresponding di- 



agonalized Hamiltonian obtained from H by a similarity 
transformation. Also z is the coordination number of the 
lattice and Tr denotes the trace. For a square (honey- 
comb) lattice, z — 4 (3). Note that the energy spectrum 
of a lattice is bounded between —z and z. By expanding 
|x| = \Hq/z\ into a Chebyshev series T2k(x), it follows 
that 

N = - + -Lip3T^(x). 

fe>l 



Also, using the identity 



T 2k (x) = (-l) k kJ2(-l) 



{2x) 2 \ 



1=0 



[k- 1)1(21)1 

we obtain 

fe>l 1=0 

where 

_ -4 , (k + l-l)\ 
1 [ z 2> (2l)l(k-l)V 

We have replaced Tr(iJQ ( ) by Tr(H 21 ), as they are equal 
to each other. 

Now let us examine the term Tr(H 21 ) more closely. 
Assuming periodic boundary conditions on the lattices, 
we have 



Tr{H 21 ) 



=< l\H 2l \l > + < 2\H 2l \2 > 



< N\H 2l \N > 



N<j\H 2l \j> 



(2) 



The Fermi-sea ground-state energy is then exactly given 
by 



£ ; T( $ ) = _£ {1 _ 2 ^_A_ [ ^r i M 2i( $)]}, (3) 



fc>i 



1=0 



where $ equals 2-7T times the flux through each lattice 
plaquette, and M 2 ;($) =< j\H 2l \j > is a moment, or 
lattice path-integral, discussed in more detail below. 

We have also studied two additional expansion 
schemes; one of them using a different set of orthogo- 
nal polynomials and the other based on a power series. 
One of them exploits the expansion of \x\ in terms of 
Legendre polynomials, P2k(x), as 



4 + E 



(2* -3)1! 



kA (2k + 2 V .l 



(4k + l)P 2k (x). 



After expressing each P2k(x) m terms of a power series, 
it is clear that 
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Et($) 

fe>i v ; i=a 

where 

-1, (2fc + 2Q! 
1 [ z 2 ' [2l)\{k-l)\{k + l)V 

The other scheme starts by rewriting Et as follows 
Et = -^Tr[(HZ) 1/2 } = -^Tr{[l + (-1 + |° )]V* }j 

where o is a constant to ensure | — 1 + Hq/o, 2 \ < 1. In 
general, a = 2\/2 for a square lattice, and a = 3/\/2 for 
a honeycomb lattice. However, for a magnetic flux near 
1/2, we find that a — 2 produces better results for both 
lattice types, since the energy density vanishes for larger 
values of E = 2\[2. By direct series expansions, we then 
obtain 

^(*) = f{-i+E A *E^p^™- (5) 

fe>l 2=0 v y 

with 

_ 1 (2fc-2)l 
fe 2 2fe -!fc [(fc- l)!] 2 ' 

Among the above three solutions for £7t($) (namely, 
Eqs. (3-5)), Eq. (3) gives the best results because of its 
relatively rapid convergence. We will return to the dis- 
cussion on their accuracy later. 

Sum over paths. — The lattice path- integral (or mo- 
ment) used in this work is 

M 2l =<j\H 2l \j>= Yl e< * e > ( 6 ) 

All 2D closed paths 

where $ c /27r is equal to the net flux enclosed by the di- 
rected closed path and 21 is referred as its order. The 
physical meaning of the above quantum mechanical ex- 
pectation value is simple. The Hamiltonian H is ap- 
plied 21 times to the initial state \j > localized at site 
j. This provides enough kinetic energy for the elec- 
tron to hop through 21 bonds, reaching the new state 
H 2l \j > located at the end of the path. The above ex- 
pectation value is non-zero only when the path ends at 
the starting site j. In our problem, quantum interfer- 
ence arises because the phase factors of different closed 
paths, or separate contributions from the same path, in- 
terfere with each other, sometimes producing cancela- 
tions in the phases. Therefore, our calculations general- 
ize the Aharonov-Bohm phase factor, obtained by having 
an electron going around a single flux-enclosing-loop, to 
the (multiply-connected) lattice case. 

The basic problem is now reduced to the computa- 
tion of the lattice path-integrals. This is a very difficult 



task, since each moment involves an enormous number of 
different loops, each one weighted by its corresponding 
phase factor. We have considerably simplified this cal- 
culation by analyzing the symmetries of the problem. In 
the next few paragraphs, we will list the most important 
symmetries involved and we will present a few examples 
of how the method works. Further details, with applica- 
tions of this method to other problems, will be presented 
elsewhere. 

We will now compute M 2 \. We consider a unit spacing 
for square as well as graphite lattices and employ the 
Landau gauge A = (0,Bx). Note that Mq equals 1 and 
is independent of the type of lattice. Also, the M 2 {s 
are gauge invariant. First let us investigate the square 
lattice. Writing the coordinates of site j as j = (m,n), 
we define an auxiliary quantity, W r (m,n), which is the 
sum over all possible paths of r steps on which an electron 
may hop from some given site to the site (m,n). From 
the definition of W r (m, n), it is evident that it obeys the 
recurrence relation 

W r+1 (m,n) = W r (m±l,n) + e* im *W r {m,n±l). (7) 

Equation (8) states that the site (m, n) can be reached by 
taking the (r + l)th step from the four nearest-neighbor 
sites. The factors in front of the W r 's account for the 
presence of the magnetic field. We can construct further 
recurrence relations successively. For example, 

W r+2 (m,n) = 4W r (m,n) + (1 + e ±4 *) 

x [e im ^W r (m ± 1, n - 1) + e' lm ^W r {m T 1, n + 1)] 

+W r (m ± 2, n) + e ±2im *W r (m, n =p 2). (8) 

Examining the action of the Hamiltonian on the state 
\j >= \m,n >, we find that 

-H\m,n>= \m±l,n> +e =F<ro *|m,n± 1 > . (9) 

Hence, by comparing Eq. (7) with Eq. (9), we obtain M 2 i 
which is just the coefficient of W r (m, n) in the recurrence 
relation for W r+2 i(m, n). This coefficient is obviously the 
sum over all possible paths which return an electron to 
its original site (m, n) after hopping 21 steps. Each path 
has a phase factor corresponding to the net flux going 
through the directed (e.g., — $ c clockwise and $ c counter- 
clockwise) path. 

It is worthwhile to notice the following symmetries 
when constructing recurrence relations and obtaining 
M 2l . 

(a) The recurrence relation for W r +i contains only 
terms W r (m + p,n± q) and W r (to — p,n± q) which sat- 
isfy the restriction : p + q = 1, 3, . . . , /, for / odd and = 
0,2, ... ,1, for I even. 

(b) The coefficients in front of the W r 's can be factored 
into two parts. Each multiplicative factor involving coor- 
dinate to is always of the form e T ' ,m3> , for W r (m+p, n±q) 
and W r (m — p,n± q). We shall refer to the rest of the 
prefactors (the part not involving to) as C r . For instance, 
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C r (m±l,n-1) = C r (mTl,n+l) = l+e*** in the above 
expression (Eq. (8)) for W r +2(m, n). These C r 's satisfy 
C r (m + p,n — q) = C r (m — p, n + g) = C r (m, + q,n — p) = 
C r (m — q, n+p) and C r (m+p, n + q) = C r (m—p, n — q) = 
C r (m + q, n+p) = C r (m — q, n — p). It can be shown that 
the latter set becomes equivalent to the former one when 
$ — > — $ (and vice- versa). 

(c) To obtain M% , it is sufficient to compute the com- 
plete recurrence relation for W r -\-i(m, n). 

We have computed the path- integrals up to M40. 
Here we list M 2 through M w : 4, 28 + 8cos$, 232 + 
144 cos $ + 24 cos 2$, 2156 + 2016 cos $ + 616 cos 2$ + 
96 cos 3$ + 16 cos 4$, 21944+26320 cos $ + 11080 cos 2$+ 
3120 cos 3$ + 840 cos 4$ + 160 cos 5$ + 40 cos 6$, for I = 
2, • • • , 10. Notice that moments with odd orders are al- 
ways zero because there is no path with an odd number 
of hops for which the electron may return to the initial 
site on a bipartite lattice. 

We now consider the hexagonal lattice, which consists 
of two interpenetrating triangular sublattices. Follow- 
ing similar techniques, we have two different formulae for 
W r +i{m, n). The proper choice between them depends 
on the sublattice to which the site (m, n) belongs. How- 
ever, both choices lead to the same results for the path- 
integrals. These have been computed up to M 60 . We list 
M 2 through Mm here: 3, 15, 87 + 6 cos $, 543 + 96 cos $, 
3543 + 1080 cos $ + 30 cos 2$, for I = 2, • • • , 10. 

In general, the path-integrals can be analytically com- 
puted by hand to any desired order through the tech- 
niques discussed above. We have computed by hand the 
moments up to the 20th order for both square and hexag- 
onal lattices. However, these and the higher order mo- 
ments can be most conveniently obtained by using com- 
puter symbolic-manipulation software. The correctness 
of the calculated moments is assured by the consistency 
of the results obtained by hand and by computer. 

Fermi-sea energy values and discussion. — After com- 
puting the lattice path- integrals, we can now proceed to 
calculate the kinetic energy of the half-filled Fermi sea. 
Recall that we have obtained the moments up to M40 
{Mqo) for the square (honeycomb) lattice. Therefore, 
by truncating the series at k = 20 (30) for a square 
(honeycomb) lattice in Eqs. (3-5), we obtain analytic 
closed-form expressions for the ground state energy as 
an explicit function of the magnetic flux. In Table I, we 
present our results for the Fermi sea energies at various 
values of the flux by using Eqs. (3,4). For the square 
(hexagonal) lattice, the results are obtained by using the 
path-integrals up to M 2 o, M 30 and M 40 (M20, M± and 
A/go) respectively. Here, instead of showing the exact nu- 
merical expressions (involving, e.g., 7r's and square roots 
of integers), we present the actual numerical values for 
easier comparison purposes. It is worthwhile to notice 
that results obtained by using moments up to M20, are 
already in excellent agreement with those obtained in 
Refs. 6,7. The values obtained by using higher order mo- 
ments are essentially identical. Also Eq. (3) and Eq. (4) 
produce almost identical results. Although the results 



obtained from Eq. (5) are not as close as those obtained 
from Eq. (3,4), they are consistent within ±0.03. 

It should be pointed out here that the tight-binding 
model H does not include the diamagnetic energy of the 
tight-binding orbitals and the reduction of the hopping 
amplitude by the magnetic field. This issue is outside the 
scope of this paper. The interested reader can find a very 
detailed analysis of these points in Ref. 12 and references 
therein. We also note that an expression for the den- 
sity of states, in terms of elliptic integrals and obtained 
through a completely different approach, has been known 
for some time. and used, for instance, in Ref. 7. Also, sev- 
eral groupstffl, including ours, have obtained results by 
other methods, including purely numerical approaches. 
Furthermore, different approaches on similar problems 
have been recently explored!!! 

In conclusion, the theory of electronic diamagnetism 
in two-dimensional lattices has been studied extensively 
due to its many applications in very diverse areas of 
physics. In particular, several computations of the Fermi- 
sea energy have recently attracted considerable atten- 
tion by many workers. We use a conceptually novel ap- 
proach that maps the problem exactly onto a Feynman- 
Vdovichenko lattice walker. More specifically, we derive 
an expression for the Fermi-sea kinetic energy at half- 
filling, as a function of a uniform perpendicular magnetic 
field, in terms of the quantum interference originating 
from the sum over 2D lattice closed paths, each loop 
weighted by the phase factor corresponding to the net 
flux enclosed. The energies obtained are essentially iden- 
tical to the ones obtained through numerical calculations. 
We have shown that lattice path-integral techniques can 
be successfully applied to this system and we expect this 
approach to be applicable to many other electronic prob- 
lems. 
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TABLE I. Fermi-sea ground-state energies at half-filling for the square and hexagonal lattices and for various values of the 
flux. The first three (4th to 6th) rows present results obtained from Eq. (3) (Eq. (4)), i.e., by using a Chebyshev (Legendre) 
series expansion. The numbers in parentheses indicate the highest order of the path-integral used. They are compared with 
the numerical results (Num.) from Refs. 6,7. 

Square lattice 



•I> 





1/8 


1/6 


1/4 


1/3 


3/8 


1/2 


C (20) 


-0.8114920 


-0.8254256 


-0.8288598 


-0.8592316 


-0.8569940 


-0.8801229 


-0.9583405 


C (30) 


-0.8109902 


-0.8231301 


-0.8347727 


-0.8599869 


-0.8576158 


-0.8775567 


-0.9581710 


C (40) 


-0.8108108 


-0.8252121 


-0.8361273 


-0.8588387 


-0.8574827 


-0.8774209 


-0.9580550 


L (20) 


-0.8114462 


-0.8250180 


-0.8286112 


-0.8603026 


-0.8567193 


-0.8799075 


-0.9584124 


L (30) 


-0.8109779 


-0.8234147 


-0.8348405 


-0.8600546 


-0.8576395 


-0.8777049 


-0.9582486 


L (40) 


-0.8108047 


-0.8254116 


-0.8362933 


-0.8587258 


-0.8575003 


-0.8774703 


-0.9580659 


Num. 


-0.811 


-0.826 


-0.835 


-0.859 


-0.857 


-0.880 


-0.958 


$ 







Hexagonal lattice 

1/4 


1/3 




1/2 


C (20) 


-0.7869927 




-0.7506642 




-0.7489389 




-0.7527022 


C (40) 


-0.7872330 




-0.7527290 




-0.7506511 




-0.7538837 


C (60) 


-0.7872775 




-0.7530423 




-0.7509664 




-0.7537718 


L (20) 


-0.7869492 




-0.7512373 




-0.7486450 




-0.7524125 


L (40) 


-0.7872308 




-0.7527266 




-0.7506855 




-0.7537969 


L (60) 


-0.7872783 




-0.7530359 




-0.7510000 




-0.7537457 


Num. 


-0.787 




-0.753 




-0.751 




-0.754 
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